Riding a wild horse: Majorana fermions interacting with solitons of fast bosonic fields. 
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I consider a class of one-dimensional models where Majorana fermions interact with bosonic 
fields. Contrary to a more familiar situation where bosonic degrees of freedom are phonons and 
as such form a slow subsystem, I consider fast bosons. Such situation exists when the bosonic 
modes appear as collective excitations of interacting electrons as, for instance, in superconductors 
or carbon nanotubes. It is shown that an entire new class of excitations emerge, namely bound 
states of solitons and Majorana fermions. The latter bound states are not topological and their 
existence and number depend on the interactions and the soliton's velocity. Intriguingly the number 
of bound states increases with the soliton's velocity. 



I. INTRODUCTION 



Models describing Majorana fermions have received a 
lot of attention, mostly in the context of quantum com- 
putation. In the process of the discussion several models 
with interesting common features have emerged. In all 
these models Majorana fermions interact with solitons of 
bosonic fields representing collective degrees of freedom 
of electronic systems (usually phase fields of supercon- 
ducting order parameters) . Although models of interact- 
ing solitons and fermions have been considered before and 
there is an extensive literature on the subject, there are 
some new features which merit a discussion. The most 
important new feature shared by the models in question 
is that Majorana fermions are slow in comparison to the 
bosonic modes. This is opposite to a more familiar situa- 
tion of the Peierls-Froelich model of polyacetylene where 
a slow optical phonons interact with fast fermions. 

As I will demonstrate, in this situation new branches of 
fermion-soliton bound states emerge some of which exist 
only in a finite region of momentum space. 



II. MODELS 
A. p-wave Josephson junction 

This model was formulated by Grosfeld and Stern [1] 
who considered a long insulating one-dimensional Joseph- 
son junction between two p-wave superconductors. As it 
is well known, p-wave superconductors have zero energy 
Majorana modes as boundary states. When the bound- 
ary is extended, as in the case of a long junction, these 
modes propagate along the edge with velocity ^ A, 
the superconducting gap in the bulk. A conventional 
long Josephson junction is described by the sine- Gordon 
model; the p-wave one acquires an additional term in 
the Hamiltonian corresponding to Majorana fermions. 
The resulting model of the junction has the following 



Lagrangian density: 
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where r, I are right- and left-moving Majorana fermion 
fields propagating along different sides of the junction 
and $-field represents the phase difference between the 
two superconductors. The parameters of the model are 
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where d is the thickness of the barrier, A^, is the London 
penetration depth, is the hight of the junction and Aj 
and TO are determined by characteristics of the junction. 
The magnetic field changes the argument of the sinus to 



$/2 $/2 - eBdx/ch. 
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Model H]) is a generalization of the Super sine-Gordon 
(SSG) model; for the latter case 



TO = v/ Xj, c = V. 
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so that the model is Lorentz invariant. The SSG model is 
integrable 0] and its excitations and S-matrix are known 
In particular, it is known that solitons of SSG model 
obey non-Abelian statistics. 

Since all candidates for p-wave superconductivity have 
small critical temperatures corresponding to small v, the 
ratio v/c is also likely to be very small (however, it may 
be increased somewhat by putting the junction on top of 
a dielectric with large e). 



B. Combination of spin-orbit interaction and 
superconductivity 

Although this model is not qualitatively different from 
the previous one, there are certain features which make 
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its material realization easier. Naniely, the chiral Majo- 
rana modes emerge here not from a p-wave supercon- 
ductor, which remains a rather exotic object, but by 
other means. Namely, Majorana fermions may emerge 
in a semiconductor with a strong spin-orbit interac- 
tion subject to external magnetic field brought into con- 
tact with an s-wave superconductor (see, for example 
p],[5]). Following 4] I write down the Hamiltonian of 
two-dimensional film of such material: 

H = j ^+{x,y)n^{x,y)dxdy, (5) 

H = [p^/2m - pl]t^ + u{py(T''- - Px(J^)tz + 

where act in particle-hole and cr" in spin space re- 
spectively. As was shown in Q, when function V{y) — 
B{y) + A(?/) changes sign Majorana zero modes emerge. 
The corresponding operators of right- and left moving 
modes are made of combinations 
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C. Half filled carbon nanotube 

In the author and Nersesyan derived an effective 
field theory for armchair carbon nanotubes using the 
bosonization approach with a partial refermionization. 
At half filling the model is similar to ([T]), but the num- 
ber of Majorana fermions is not one, but 6 with different 
mass parameters TOq such that 



2toc + 7111 + "^0 = 0- 
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The role of the bosonic field $ is played by the total 
charge field <i>c. The smallness of parameter /3 and a 
large value of c/v originate from the unscreened Coulomb 
interaction. The Lagrangian density is 



^ X! Xaijodt + v-fidj^'jxa-V, 

a=-2 
V = 

icos[V47r$c] 
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These modes are spatially separated: the mode r emerges 
at the edge with dV/dy > and Z-mode at the edge 
dV/dy < 0. These edges are boundaries between the 
topological insulator and the superconducting state. Be- 
ing projected onto these modes the Hamiltonian density 
becomes 



Ho 



dx{—rdxr + Idxl)- 



(8) 



Since the spin-orbit interaction is typically much greater 
in magnitude than the p-wave order parameter, one can 
increase the ratio u/c. In Q the authors cite w = 7.6 x 
lO^cm/sec for InAr. With Al - lO^A and d - lOA one 
can get u/c ~ 10~^. 

Now following Q consider a situation when a narrow 
superconducting region is sandwiched between two topo- 
logical insulators. Then instead of one bosonic mode as 
in ID) we will have more. Namely, if the superconducting 
region between two topological insulators is sufficiently 
narrow, we have the following action 



-J COs(0a - (l>b) + 
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where (pa.b are superconducting phases on the left and 
right from the superconducting strip and (jjm is a phase 
on the strip. Therefore there are two independent bosonic 
modes. 



where 70,71 are Dirac gamma matrices, K << 1 is 
the Luttinger parameter determined by the long range 
Coulomb interaction, $c is the total charge field and Xa 
are Majorana fermions made of chiral components of 
(a = -2, -1) and fields (a = 0, 1, 2, 3) 0. The V- 

term represents the leading interaction generated by the 
Umklapp processes; the interactions between the Majo- 
rana fermions are small in comparison. The symmetry 
of mi) is U(l)xU(l)xSU(2)xZ2. The Majorana modes 
with a ~ 1,2,3 realize the S—1 representation of the 
SU(2) group. 

For X << 1 $c is essentially a classical field and its 
dynamics is determined by the equation 
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At if < < 1 one can neglect coordinate dependence of the 
mass term and treat it as a constant for the purposes of 
calculation of the fermion average in ()13p . The result is 
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where $ = ^/IGn^c and is calculated with the loga- 
rithmic accuracy. This description is valid for excitations 
moving with velocities < v. So we see that the solitons 
of <I>c are solutions of the sine-Gordon equation. As far 
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as the fermionic excitations are concerned, they are de- 
termined by the same equations as for model ([1]). In the 
notations of ^ we have 



1/A,; 
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As is shown the subsequent Section, the fermions have 
bound states with the sohtons such that the number of 
finite energy bound states in each channel is 



No,a = rriaXj/v 
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It is clear that for the Lorentz invariant case K = 1 there 
are no finite energy bound states. More than that, they 
appear only if the Coulomb interaction is quite strong. 



III. SEMICLASSICAL ANALYSIS 

In all models described above the bosonic action is of 
the sine-Gordon type. The sine-Gordon subsystem has 
two types of excitations; kinks and breathers. Kinks 
strongly interact with the Majorana fermions since the 
latter ones create bound states with kinks. There are 
two types of bound states: one type is the Majorana 
zero modes which modify the kink's quantum numbers 
and the other are massive ones. Below I do the analy- 
sis for model ([TJ, generalizations for models (|9llip are 
straightforward. In particular, in model (jlip c — v / K. 

It is convenient to introduce new fermionic fields: 



Xi-(r + 0/V2, X2 = {r-l)/V2 



(17) 



I consider a static kink configuration first. Then for a 
kink centered around xq the fermion operators can be 
represented by the mode expansion (IT51) [8i]. The eigen- 
functions satisfy 



= i-v^dl + W'^ - vd^W)^IJi (18) 



where W = mcos($/2). 







In the static case u = this gives rise to the following 
potential: 



W = mtanh(a;/Aj). 



(21) 



Substituting it in and using the results from 0, I 
obtain the following eigenvalues for the bound states: 

n = 0,1,... < No mXj/v. (22) 

and the eigenfunctions are 

V'l = (coshO^"^''+"^ X (23) 
F(^~n,2No-n + l,No-n + l;-^^^y 

Notice that in the Lorentz invariant case ^ No = 1 and 
the only bound state is the topological one n ~ 0. 

Now let us consider the case of moving soliton u ^ 0. 
Let us introduce new coordinates: 

x' = j{x — ut), t' = j{t — xu/v^), 
-1^[l~{u/vf]-^''^. (24) 

This Lorentz transformation leaves the fermionic action 
invariant and puts us in the reference frame of the moving 
soliton. The mass term in becomes 
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Notice that in the Lorentz invariant case c = v the scale 
does not change. However, if t; < c, the soliton size 
increases and, as a consequence, the number of bound 
states N also increases: 
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This is a somewhat unexpected result. The energy in this 
reference frame is given by (|22p with A^o replaced by N. 

In the laboratory reference frame the energy and mo- 
mentum of the fermionic part of the bound state are (I 
set c oo): 
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A general single-soliton solution is 



To obtain the total energy and momentum of the kink 
and the bound state one has to add the energy and mo- 
mentum of the kink: 
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FIG. 2: The spectrum e = No[E„ip) - Mk]/2m ^ for No = 
5 and n — 1, 2, 3. The fermion velocity is set v = 1. 
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FIG. 1: The spectrum e = No[E{p)-Mk]/2m ^ for A^o = 5 
and n — 6. The spectrum exists only for momenta larger than 
critical as explained in the text. The fermion velocity is set 
v = l. 



Consider n = mode first. For zero velocity this mode 
always exists, but for finite velocities its existence is re- 
stricted by the condition u < v. Assuming that c/v is 
so large that the kink's momentum is much smaller than 
the fermionic part and its dispersion is slow, we can set 
Ek — Mk, Pk = 0. Then from ([^5]) we extract its disper- 
sion: 

Eip) = ^Ml + {cp)\ (29) 
c 

Thus the zero energy bound state always exists, though 
in a limited region of the Brillouin zone. 

Now consider the finite energy bound states. They are 
not topological and their existence is conditional. Some- 
what unexpectedly (l?7|) shows that the conditions for 
their existence improve when the kink's velocity increases 
towards v. There are solutions with iVo[l — > 
n > Nq existing only for finite velocities (momenta). It 
is illustrated on Fig. 1. 

To derive the dispersion of the bound states we have 
to take into account the fact that their energies and mo- 
menta are sums of ([77]) and If we assume that 

2ml No » Mk{v/cf, (30) 

then the inertia of the kinks is very small and their contri- 
bution to the total momentum can be neglected in com- 
parison with the momentum of the bound state (j27p . As 
a result one gets the picture of the dispersion depicted 
on Fig. 2. 



IV. QUANTUM NUMBERS AND 
CORRELATION FUNCTIONS 

As we see from ([T9|) the operators of fermion-kink zero 
modes 7o(a::o) compose a Clifford algebra. For models 
with several species of Majorana fermions, such as (fTT|) . 
7q create a spinor representation of the corresponding 
group (for (HII) the group is 0(6) SU(4) ) and the 
bound states of solitons and Majorana fermions trans- 
form according to this spinor representation. For the 
case of half filled carbon nanotube these excitations carry 
the same quantum numbers as the original fermions and 
therefore are quasiparticles. The situation with n ^ Q 
bound states is quite different. They transform according 
to the vector representation of the corresponding group 
and therefore can be created only by pairs of fermionic 
operators. 

In order to get a better grip of the picture, let us 
consider the model of carbon nanotube. The symmetry 
group of model (HI]) is U(l)xU(l)xO(3)xZ2. As an ex- 
ample of a local field having nonzero matrix elements be- 
tween the vacuum and the aforementioned bound states 
we have 

gix/4¥v'c cos[VA^ipf] - Rl^Rl^ + Rt^Rti. (31) 

where Lp^.f are right- moving components of the corre- 
sponding bosonic fields, c labels the total charge, / la- 
bels the relative one and 1, 2 label positions of two Dirac 
points in the Brillouin zone of carbon nanotube. The first 
exponent creates two right-moving solitons in the charge 
sector. One soliton creates a bound state with a Majo- 
rana fermion from /-sector (asymmetric charge) and the 
other soliton remains unbounded. From this example one 
can see that the bound states can be observed only as 
parts of continua. In the above example the continuum 
consists of a " naked" soliton and a soliton- fermion bound 
state. Existence of "naked" solitons, i.e. ones which do 
not carry any fermionic modes is guaranteed by the fact 
that c > V and so their is plenty of room in momentum 
space for solitons which velocity exceeds the one of the 
fermions and those, as we know, do not create bound 
states. 



CONCLUSIONS a limited region of momentum space. 



This paper demonstrates that in field theories with- 
out Lorcntz invariancc (quite a common thing in con- 
densed matter physics) one has to expect appearance of 
new types of bound states, some of them existing only in 
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